Multi-strategy synthetized equilibrium optimizer and application

Background Improvement on the updating equation of an algorithm is among the most improving techniques. Due to the lack of search ability, high computational complexity and poor operability of equilibrium optimizer (EO) in solving complex optimization problems, an improved EO is proposed in this article, namely the multi-strategy on updating synthetized EO (MS-EO). Method Firstly, a simplified updating strategy is adopted in EO to improve operability and reduce computational complexity. Secondly, an information sharing strategy updates the concentrations in the early iterative stage using a dynamic tuning strategy in the simplified EO to form a simplified sharing EO (SS-EO) and enhance the exploration ability. Thirdly, a migration strategy and a golden section strategy are used for a golden particle updating to construct a Golden SS-EO (GS-EO) and improve the search ability. Finally, an elite learning strategy is implemented for the worst particle updating in the late stage to form MS-EO and strengthen the exploitation ability. The strategies are embedded into EO to balance between exploration and exploitation by giving full play to their respective advantages. Result and Finding Experimental results on the complex functions from CEC2013 and CEC2017 test sets demonstrate that MS-EO outperforms EO and quite a few state-of-the-art algorithms in search ability, running speed and operability. The experimental results of feature selection on several datasets show that MS-EO also provides more advantages.


INTRODUCTION
Traditional mathematical methods have strict requirements for optimization problems (OPs) (Hashim et al., 2019).They cannot get good results in solving most OPs.Inspired by natural phenomena and biological behavior (Ghasemi et al., 2023a), researchers have proposed many meta-heuristic algorithms (MAs) to solve OPs better.They include genetic algorithm (GA) (Zhou et al., 2021), simulated annealing (SA) (Kirkpatrick, Gelatt & Vecchi, 1983), particle swarm optimization (PSO) (Chen & Lin, 2009), differential evolution (DE) (2020) applied an adaptive EO (A-EO) to multilevel thresholding image segmentation.A-EO solved the problem that search agents were randomly scattered to nonperformer search agents in EO.Biller et al. (2016) proposed an improved EO through linear classification reduction diversity technique and local minima elimination method, while the variant of EO proposed by Ghasemi et al. (2023a), Ghasemi et al. (2023b) and Ghasemi et al. (2023c) used different probabilities to select equilibrium candidate solutions.They were applied to photovoltaic parameter estimation, and effectively improved the optimization accuracy and reliability.These EO variants have improved the optimization performance of EO to different degrees.Although they all verify the performances of the classical functions, there are no reports on verifying the variants of EO on CEC2013 and CEC2017 test sets.
As mentioned above, modification and addition have their own shortcomings and advantages.A hybridization of the two can compensate for their respective shortcomings, and an ingenious combination of several strategies can allow individuals to work closely together to balance exploration and exploitation.Thus, the main purpose of this article is to propose many strategies based on updating equations and discuss how they can be integrated to improve the overall performance of EO and solve some of its drawbacks.Inspired by the above description, this article proposes a multi-strategy synthetic EO (MS-EO).The graphical abstract of this article is shown in Fig. 1.The contributions of this article are as follows: • A simplified updating strategy is applied to EO to form a simplified EO (SEO) by simplifying the original updating equation to improve the operability and reduce the computational complexity.
• An information-sharing strategy is added as an updating way to SEO and two updating ways alternate in different iterations to improve the global search ability.
• A migration strategy is used to update a golden particle to get a stronger search ability.
• An elite learning strategy is implemented for the worst particle updating in the late iteration stage to enhance the local search ability.
The organization of this article is structured in the following order: the related work, the proposed algorithm (MS-EO) , the experimental results and analysis, the application of MS-EO to feature selection, conclusion and future work.

Related work
EO is inspired by the law of dynamic mass balance.Its mathematical updating model is as follows: where C i represents the concentration of the i-th particle.C e is the concentration of the elite particle selected randomly from the elite group.F is an exponential term, G represents the generation rate, and λ is the turnover rate and is often taken as a random vector with each component distributed in [0,1].V is the volume and is usually regarded as a unit.+ means term by term multiplication.i ∈[1, N ], and N is the population size.
In order to maintain the diversity of the population and the direction of evolution, in the iterative process, some of the best individuals in the current population will form the Experiment results on CEC2013 and CEC2017 test sets demonstrate that MS-EO outperforms EO and quite a few state-of-the-art algorithms.

Test on 58 complex functions
Experimental results of feature selection on 12 datasets show that MS-EO also provides more advantages.
Application to feature selection elite group to pass to the next generation.In this article, the elite group is composed of five elite particles, which are as follows: The first four are the best-so-far particles' concentrations in the current population, the first C α , the second C β , the third C δ and the fourth C γ best one, respectively.The last particle C ϕ is the arithmetic mean of the four best particles' concentrations.
The exponential term F is expressed as: where a 1 is a tuning parameter, sign is a sign function, and r is a random vector in [0, 1] like λ.ε is a dynamic tuning parameter and its expression is as follows: where t is the current iteration number, T represents the maximum iteration number, and a 2 is also a parameter.F is also a dynamic random parameter vector and its component value is shown in Fig. 2. Fig. 2 only illustrates the value changes on the six components of F with T equal to 1,000.From Fig. 2, the red line represents the value change of F with Generation rate G can be calculated as: where R G is the generation rate control parameter and is determined by P G .P G is Generation Probability and a parameter.r 1 and r 2 are random numbers in [0, 1].From Eqs. ( 5), ( 6) and ( 7), G is a zero vector when r 2 is less than P G .From Eq. ( 1), the first term on the right of the equation is the concentration of the elite particle, and the second term is the difference between the elite particle's and the i-th particle's concentration.The third term is the generation item.Its pseudo-code is shown in Algorithm1 EO.From the above description, EO has the following good features.(1) Unique search mechanism: Each particle's concentration is updated in a unique updating way including the elite particle selected randomly from an elite group.(2) High degree of information guidance: The three items on the right of Eq. (1) all contain elites guiding each particle.(3) Strong local search ability: The guidance of the elite particles makes EO evolve in a good direction throughout the iteration process.Each component value of F is small in the late iteration stage.The difference between the elite particle and other particles is small, and the value of G is small.Meanwhile, the second term on the right side of Eq. ( 1) takes up less weight.Moreover, the particle is closer to the elite particle when R G is 0. These help EO search more accurately.(4) Global search ability.Each component value of F and the difference between the elite particles and other particles are large in the early search stage.Each component value of G is also large when λ randomly takes a small value and R G is not 0. At this time, the value of the third term is also larger, which helps the current particles search the entire search space.
However, there are some shortcomings when EO is used to solve complex OPs.(1) From Eq. ( 1), the particles in EO learn from the elite particles, and the generation term contributes to the accurate search.This will make EO easy to fall into local optima and have poor exploration ability.(2) From Eq. ( 1), the particles in the population only communicate with elite particles.This results in a low degree of information sharing and poor diversity in the population.Moreover, EO only learns from elite particles and is also monotonous, which has no other way of guiding and is not conducive to the further evolution of the population.It leads to insufficient search ability.(3) The computational load of the exponential function is higher.Moreover, λ is a random vector.It calculates the exponential term F more complicated.(4) EO has many parameters, such as λ, F , r 1 , r 2 , r, V, a 1 , a 2 and P G .Although the parameters increase the flexibility of EO, it leads to poor operability of EO.

MATERIALS & METHODS
As mentioned in the previous section, EO has a number of advantages when it comes to solving classic OPs (Faramarzi et al., 2020).However, there are still drawbacks when solving complex OPs.Therefore, this article proposes a Multi-strategy Synthetized EO (MS-EO) to cope with the above problems of EO.

Simplified updating strategy
In order to reduce the computational complexity and improve the operability of EO, a simplified updating strategy is used in EO to form a simplified EO (SEO).Its expression is as follows: Figure 3 shows the values of simplified F .From Eqs. (8), ( 9) and (10) and compared with EO, the simplified updating strategy has the following differences: (1) The simplified updating strategy removes the parameters λ.Each component value of F in Fig. 3 is more regular and the curves are smoother compared with Fig. 2. Therefore, the simplified updating strategy reduces the computational complexity of EO.Full-size DOI: 10.7717/peerjcs.1760/fig- 3 (2) The generation term in Eq. ( 8) does not have a value of 0 compared with Eq. ( 1).The computational load descends to a certain extent.
(3) From Eq. ( 8), the simplified updating equation highlights the elite particles even more and also highlights the exploitation.
In summary, the simplified updating strategy is the modification of the updating equation.It can effectively reduce the computational complexity of EO and improve its operability.However, the simplified updating strategy does not change the overall search pattern of EO.Therefore, it still has insufficient search ability.

Information sharing strategy
Whether the information of individuals in the population is fully utilized or not is crucial to the search ability of a population-based MA.The performance of the MA will be significantly improved if the information can be fully utilized in the search process (Zhang & Yang, 2021).In order to enhance the information sharing and the global search ability of EO, this article introduces an information sharing strategy.The mathematical model of the information sharing strategy is as follows: where C a and C b are the concentration of two particles selected randomly from the current population, i =a =b and inspired by Zhang et al. (2020b), the scale factor f r adopts the dynamically adjusted sinusoidal approach.From Eqs. ( 11) and ( 12), a new solution is generated based on the combined action of three different solutions.η is the difference between two random solutions, which realizes the information sharing among particles in the population, because all the particles in the population may be selected with the increase of the iteration number t, it makes full use of the information of the whole population and helps increase the population diversity.
The simplified updating strategy still maintains the strong exploitation ability of EO.In order to enhance the global search ability of EO, the information sharing strategy and the simplified updating strategy are combined by a dynamical adjustment strategy to form a simplified sharing EO (SS-EO).The flow chart is shown in Fig. 4. According to Eq. ( 10), ε is a dynamically adjusted parameter of the simplified strategy that linearly decreases from 1 to 0.
From Fig. 4, the current particle is updated with the simplified updating equation when r 2 is greater than or equal to ε.Otherwise, the current particle uses the information sharing updating equation to update.In the early iteration stage, the value of ε is larger.Most of the particles are updated by the information sharing strategy, and a few of them are updated in a simplified way.The value of ε is small in the late iteration stage.Most of the particles are updated by the simplified updating way with strong exploitation ability, and a small number of particles are updated by the information sharing strategy.

Golden particle migration strategy
The golden section is a proportional relationship in mathematics.It enables us to reasonably arrange fewer test times.The migration operator in BBO (Zhang et al., 2019) is to share the characteristics of good individuals with poor individuals.It can effectively improve the search ability of poor individuals, so as to enhance the search ability of the whole population.Inspired by the above and in order to stimulate and enhance the search ability, a golden particle migration strategy is introduced SS-EO to form a golden SS-EO (GS-EO).The mathematical model of the migration strategy is as follows: where C represents the j-th value of the concentration of the golden particle, C is the j-th value of the example particle, q ∈[1, g -1], and j ∈[1, D ]. g is the index of the sorted particles from the best to the worst according to the fitness values of all the particles.ceil is round-up function.The illustration of the golden particle migration strategy is shown in Fig. 5. From Eqs. ( 14), ( 15) and ( 16), an example pool is composed of the particles that are better than the golden particle.From Fig. 5, each dimension value of the golden particle comes from the corresponding dimension value of the example particle selected randomly from the example pool.The golden particle migration strategy is individual-based addition.The golden particle requires D example particles.At the same time, each dimension of the golden particle is derived from its example particle also makes it evolve in a good direction.The golden particle gradually becomes excellent as the iteration number increases.And it may be selected as the learning object during the iterative process by the simplified updating strategy and information sharing strategy to participate in the updating for other particles.This further stimulates the search for other particles, so as to improve the search ability of the algorithm.

Elite learning strategy
Educators often raise the grades of poor students by letting excellent students help poor students, so as to improve the overall grades of the class.Inspired by this and in order to enhance the local search ability, an elite learning strategy is used for the worst particle in the late iteration stage to form MS-EO.The mathematical expression of the elite learning strategy is as follows: where C w is the concentration of the worst particle in the current population.From Eqs. ( 17) and ( 18), the worst particle learns from the elite particle to make up for its shortcomings.And it is used in the late stage.Therefore the elite learning strategy belongs to individual-based addition and iteration-based addition.In the late iteration stage, the component value of θ is smaller and the value of f r fluctuates widely.When the value of f r is small, it is helpful to the local search of the algorithm.It can avoid the algorithm falling into local optima when the value of f r is large.Moreover, the worst particles learn from the elite particles.It greatly enhances the exploitation ability of the worst particle, so as to enhance the search ability of the whole population.

Framework of the proposed algorithm
MS-EO skillfully synthesizes many strategies, such as simplifying, information sharing, golden particle migration and elite learning, and embeds them into EO according to their respective characteristics.The flow chart of MS-EO is shown in Fig. 6.
From Fig. 6, the information sharing strategy is mainly used to update the most particles' concentrations in the early iteration stage.In the late iteration stage, most particles are mainly updated in a simplified way.For the golden particle, the migration strategy is adopted to make it originate from the example particles.For the worst particle, the elite learning strategy is adopted in the late iteration stage to further improve the optimization level of the whole population.

Complexity analysis
On time complexity, it is assumed that the updating equation of EO requires L 0 time unit to generate a variable of a new solution, while the updating equations of the simplified updating strategy, information sharing strategy, golden particle migration strategy and elite learning strategy require L 1 , L 2 , L 3 , and L 4 to generate a variable of a new solution, respectively.According to the described above, the simplified updating strategy removes some components in EO.Eq. ( 8) has fewer operations compared with Eq. (1).Thus, Eq. ( 8) takes much less time to generate a new solution than Eq. ( 1).The updating equations of the remaining three strategies are all simpler than Eq. ( 1).Therefore, L 1 , L 2 , L 3 , and L 4 are much less than L 0 .MS-EO has lower time complexity.

Start Set the parameters T, N and initialize the population randomly
Evaluate the fitness value of each particle, select C α , C β , C δ and C γ , calculate C φ , construct an elite group and determine the golden section particle index Limit the boundaries of each new solution, evaluate the fitness values in parallel, and update the population using greedy selection and the elite group Elite learning Eq. ( 19) Information sharing Eq. ( 13)

Simplified EO
Eq. ( 9) Golden particle migration Eq. ( 15 For space complexity, the space to be allocated for EO includes current population, elite pool and fitness value.MS-EO does not allocate additional storage space.Its space complexity is the same as that of EO.

Experimental environment and evaluation standard
In order to verify MS-EO, a large number of experiments are performed on the complex functions from CEC2013 and CEC2017 test sets.Their details can be found in Liang et al. (2012); Awad et al. (2016).CEC2013 test set includes 28 functions, among which f 1 ∼f 5 are unimodal functions, f 6 ∼f 20 are basic multimodal functions, and f 21 ∼f 28 are composition functions.CEC2017 test set contains four types of functions: unimodal functions (f 1 ∼f 3 ), simple multimodal functions (f 4 ∼f 10 ), hybrid functions (f 11 ∼f 20 ) and composition functions (f 21 ∼f 30 ).Compared with CEC2013, CEC2017 is more complex, such as more function types and more shifted functions.All experiments were conducted in MATLAB (R2014a) environment running on a PC with 3.4 GHz i7-3770 CPU and 8GB RAM.The statistical analysis software was IBM SPSS Statistics 19.
This article adopts statistics such as mean value (Mean), standard deviation (Std), ranking (Rank), and average ranking (Ave.Rank) to evaluate the optimization performance of an algorithm (Zhang et al., 2019).The algorithm with a smaller Mean value has better performance for a minimum problem in this article.The ranking rule is as follows, On each function, the less the Mean value obtained by the algorithm is, the higher the ranking is; if the Mean values by the algorithms are the same, the higher the ranking is and the less the Std value is; if the Mean and Std values of the algorithms are the same, then the algorithms rank the same.The best data in the table is bold.
For a fair comparison, the common parameter settings of all the algorithms are the same, such as the maximum number of function evaluations (M nfe ), the number of independent runs (N run ), and D. According to the best recommendation of (Liang et al., 2012;Awad et al., 2016), M nfe is set to D ×10000 and N run is 51.Different parameter settings of an algorithm will result in different experimental results, so other parameter settings of the comparison algorithms are set by the best settings from the corresponding references shown in Table 1.In addition, all algorithms assume that the experiment is carried out in an environment without interference from other factors.

Comparison with its incomplete algorithms
In order to demonstrate the contribution of each new strategy of MS-EO, this subsection compares MS-EO with EO, SS-EO and GS-EO on the 30-dimensional functions from CEC2013 test set.In order to explain the problem briefly, this subsection selects the experimental results of three representative functions to illustrate.The three representative functions are unimodal function (f 3 ), basic multimodal function (f 7 ) and composition function (f 22 ).Fig. 7 provides a performance graph of MS-EO and incomplete variants on the 3 representative functions and an average ranking chart on all the functions,   respectively.The parameters settings for MS-EO and the incomplete variants are set to the same.From Fig. 7D, the corresponding ranking order is MS-EO, GS-EO, SS-EO and EO on all the functions.From Figs. 7A, 7B and 7C, the information sharing strategy, the simplified update strategy and the gold particle migration strategy all improve the performance of EO to a certain extent, and the adoption of the elite learning strategy for the worst particles in the late iteration stage can effectively enhance the local search ability of SS-EO, which indicates that the optimization performance of MS-EO is better than that of GS-EO, especially on f 7 .It shows that every improvement of MS-EO is effective.
In Fig. 8, the height of each color block represents the times obtained by an algorithm in ranking the x-th.For example, ' 1 ' represents ranking the first.The higher the height of the color block is the more times obtained by the algorithm in ranking first.It is obvious that these three color blocks (' 1 ', ' 2 ' and ' 3 ') have the highest height in MS-EO.Therefore, MS-EO has the best optimization performance.From Table 2, on Ave.Rank, the value of EO is 6.04, ranking the eighth among all the 10 algorithms.MS-EO ranks first with Ave.Rank of 3.29.It indicates that the improvements proposed in this article greatly enhance the optimization performance of EO.

Performance comparison on the CEC2017 test set
In order to further verify the ability of MS-EO to solve more complex problems, many experiments are conducted on the more complex functions from the CEC2017 test set.Some representative algorithms are selected for comparison.These algorithms include EO, DQL-SFLA, Sa-DE, ME-GWO, A-EO, HFPSO, HBBOG, DS-PSO and MSS-CS.The parameters of the algorithms are shown in Table 1.Tables 3 and 4 provide the comparison results of MS-EO and these representative algorithms on the 30-and 50-dimensional functions from CEC2017 test set, respectively.Among them, the data of Sa-DE is directly taken from Tang (2019).
From Table 3, on Ave.Rank, EO is 7.33 and ranks ninth.It is better than HFPSO (8.60).On unimodal functions, MSS-CS ranks first on f 1 , MS-EO ranks the first on f 2 and f 3 , and other comparison algorithms obtain no times ranking the first.It means that MS-EO has strong exploitation ability.DS-PSO obtains four times ranking first on the multimodal functions.It shows the strong exploration ability of DS-PSO.However, MS-EO can obtain one time ranking first and four times ranking second on the multimodal functions.It indicates that MS-EO also has strong exploration ability compared with the other comparison algorithms.MS-EO won six times and three times ranking first on the hybrid and composition functions, respectively.However, from the overall ranking,      The results on the 50-dimensional functions are similar to those on the 30-dimensional functions.From Table 4, EO still ranks ninth and is only better than HFPSO.On the unimodal functions, MS-EO, Sa-DE and MSS-CS obtain one time ranking first, respectively.DS-PSO shows strong optimization ability on the multimodal and composition functions.However, MS-EO can obtain 6 times ranking first, and 13 times ranking the second.Moreover, the average ranking of MS-EO is 3.37 and it ranks first.It indicates that MS-EO also has the best optimization performance on the 50-dimensional function.Therefore, MS-EO shows stronger competitiveness on CEC2017 test set.In conclusion, the synthetization of many search strategies into EO is effective and MS-EO has stronger scalability and universality.

Statistical analysis
In order to verify MS-EO, this subsection performs Wilcoxon signed rank test with Bonferroni-Holm (Zhang et al., 2020a;Zhang et al., 2020b) correction based on the complex functions (D = 30 and D = 50) from CEC2017 test set (60 cases in sum).The purpose is to detect whether there is a significant difference between the two samples.Let H0 be the null hypothesis, that is, there is no significant difference between the two samples.The significance level (b) is set to 0.05 based on the best recommendation (Derrac et al., 2011).Table 5 provides Wilcoxon signed rank test results based on Table 3 and 4. R + represents the positive rank, which is the total rank of the problems in which MS-EO is better than the comparison algorithm.R − is the negative rank, which is the total rank of the problems in which MS-EO is inferior to the comparison algorithm.The corresponding rank is equally divided into R + and R − when the optimization performance of the two algorithms is equal.The p-value is calculated from R + and R − .The Bonferroni-Holm method adjusts b in descending order of p-value.H 1 , H 2 , . . ., H k -1 are corresponding hypotheses.The hypothesis is rejected when b/u < p-value, otherwise it is accepted.n/w/t /l respectively represent the number of optimization functions is n and the number of functions that MS-EO is better than, equal to and inferior to the comparison algorithm.
From Table 5, the p-values of MS-EO and the comparison algorithms are all less than b/u, and the corresponding hypothesis is rejected.It shows that there are significant differences between MS-EO and the comparison algorithms.In particular, the function number that MS-EO is better than EO is 54.Moreover, the function number that MS-EO outperforms A-EO is 52.And MS-EO is significantly superior to any of the comparison algorithms on CEC2017 test set.a fixed time if P i of the algorithm is higher, and its performance is better.The mathematical model of P i can be formulated by Eq. ( 19): where n ∈[1, m ], and m is the sum of the number of problems.On the n-th problem, M E and M S are the minimum values of the average error value and average time of all algorithms, respectively.E m and S m the average error value and average time obtained by an algorithm on the n-th problem, respectively.c 1 is the average error weight and c 2 is the average time weight (c 1 + c 2 = 1 and c 1 , c 2 ∈[0, 1]). Figure 9 shows the P i values of MS-EO and the contrast algorithms under the conditions of 30-dimensional and 50-dimensional functions in the CEC2017 test set.In addition, m = 60, the weight (w) is set to the numbers between 0 and 1 in steps of 0.2 according to the best recommendation from Gupta et al. (2020), c 1 = w and c 2 = 1-w.From Fig. 9, the value of P i of MS-EO is greater than the comparison algorithms at any time.Moreover, the value of P i of MS-EO is greater than that of EO and A-EO regardless of on CEC2017 test set.It shows that MS-EO can get more accurate solutions in a fixed time on CEC2017 test set.This further indicates that the performance of MS-EO is more stable.

Convergence analysis
To validate the convergence quality of MS-EO, this subsection performs some convergence tests on the 30-dimensional functions from the CEC2017 test set.Four representative functions are selected from the CEC2017 test set to illustrate the problem briefly, namely, unimodal function f 3 , multimodal function f 9 , hybrid function f 12 and composition function f 30 . Figure 10  algorithms.From Fig. 10, the convergence speed of MS-EO is much faster than that of the comparison algorithms on f 3 , f 9 and f 12 .On f 30 , the convergence of MS-EO is not as obvious as on the first three functions.However, MS-EO converges faster than any of the comparison algorithms as the iteration number increases, and can find a more accurate solution.On f 3 and f 12 , the convergence curve of DQL-SFLA drops suddenly and converges to a smaller solution in the late iteration stage.It is due to the fact that DQL-SFLA incorporates the quasi-Newton local search in the late stage, so that it has a faster convergence speed.In addition, the convergence speed of MS-EO is dramatically faster than that of EO and A-EO on all four functions.MS-EO has better convergence performance compared with the comparison algorithms on the 30-dimensional functions from the CEC2017 test set.

Application to feature selection
In order to verify the ability of MS-EO to solve practical problems, this section discusses MS-EO's application to feature selection.Feature selection is to find the optimal feature subset by eliminating irrelevant or redundant features from the original feature set (Hu et al., 2021).It can be defined as a discrete OP because the solution of feature selection belongs to binary.Suppose a data set has N samples and H features.B is the original feature set.Then feature selection can be regarded as finding h (h < H ) features from B while minimizing its fitness function to maximize the classification accuracy of a given classifier.Its search individual X can be expressed as: where 1 represents the corresponding feature is selected.Otherwise, it is not selected.The objective function can be defined as: where µ 1 and µ 2 are two weighting coefficients.µ 1 ∈[0,1] and µ 2 = 1 − µ 1 .According to the best recommendation (Hu et al., 2021), µ 1 is set to 0.99, and the corresponding µ 2 is 0.01.e r represents the error rate and it is calculated by K -Nearest Neighbor (KNN) classifier.KNN divides the dataset into training, validation, and testing of equal size to cross-verify each dataset (Tu, Chen & Liu, 2019).K = 5 in the experiment.Twelve datasets from UCI machine learning repository (http://archive.ics.uci.edu/ml/index.php)we selected to test the effectiveness of MS-EO on feature selection.Table 6 provides a brief description of these datasets.k, v and z represent the number of samples, features and classifications of each dataset, respectively.Among them, the distribution of D 7 (Segmentation dataset) under different features is shown in Fig. 11.From Fig. 11, four plane graphs are made by selecting two of the 19 features of D 7 .On Features 1 and 2, these samples with seven types fill the entire graph.In the remaining three feature plane graphs, the sample distribution is relatively concentrated because some sample points overlap.
Table 7 shows the comparison results of MS-EO and the comparison algorithms on the average fitness value of objective function (mean), classification accuracy (accuracy) and average number of selected features (AND).The common parameters are set as follows: For fair comparison, M nfe = 1000 for all the comparison algorithms.For MS-EO, N = 20, T = 50; For ME-GWO, N = 5, T = 100; For SR-PSO, Sin-DE and WOASA, N = 10, T = 100.The experimental results of SR-PSO, Sin-DE, WOASA and ME-GWO are all taken from Tu, Chen & Liu (2019).
From Table 7, MS-EO outperforms SR-PSO, Sin-DE and WOASA on all the datasets in terms of mean fitness.For ME-GWO, MS-EO outperforms it on all the remaining datasets except D 4 .The classification accuracy is similar to mean.On accuracy, MS-EO   8. From Table 8, the ranking order of the five algorithms is MS-EO, ME-GWO, WOASA, Sin-DE and SR-PSO according to the comprehensive ranking results in Table 8.It further proves that MS-EO can solve the problem of feature selection more effectively.

CONCLUSIONS
In order to cope with some drawbacks of EO in solving complex OPs, this article proposes a multi-strategy synthetized EO (MS-EO).Firstly, a simplified updating strategy is adopted by simplifying the original updating equation of EO to reduce the computational load and enhance the operability of EO.Secondly, an information sharing strategy based on iteration is added into the simplified EO.Most particles are updated by the strategy in the early search stage to enhance the global search ability while updated by the simplified EO in the late search stage ensuring the exploitation ability of EO.Then, a migration strategy is used  for a golden particle to update and to enhance the search ability.Finally, an elite learning strategy is used for the worst particles in the late search stage to enhance the local search ability.Experimental results on CEC2013 and CEC2017 test sets demonstrate that MS-EO is better than EO and any of some state-of-the-art algorithms.MS-EO outperforms EO on 54 of 60 functions from CEC2017 test set (D = 30 and D = 50).The running time of MS-EO accounts for 36.22% and 45.45% of EO on the 30-and 50-dimensional functions, respectively.The value of P i of MS-EO is also the highest under any condition.On feature selection, MS-EO has also obtained more convincing results.These results on both complex functions and feature selection show many updating strategies synthesized into EO are effective and expect to be extended to improvement on other MAs.In the future works, we will propose new search strategies to further improve the multi-strategy theory for MS-EO and other MAs.The improved MAs will be dealt with more practical problems.

Figure 1
Figure 1 The structure diagram summarizes four strategies of MS-EO (simplified updating strategy, information sharing strategy, golden particle migrate strategy and elite learning strategy), experimental setting and application.Full-size DOI: 10.7717/peerjcs.1760/fig-1

Figure 2
Figure 2 Values of F with t.

Figure 3
Figure 3 Values of simplified F with t.

Table 2 Friedman test results on the CEC2013 test set (D = 30).
on Ave.Rank, the value of MS-EO is 2.77.Therefore, compared with these comparison algorithms, MS-EO has strong optimization ability on the 30-dimensional functions from the CEC2017 test set.